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FOREWORD 


This  reporC  is  an  adaptation  of  Manaf  Hassan  Ali's  Ph.D.  dissertation, 
"Spectrophone  Measurements  in  Sulfur  Hexafluoride,"  which  was  accepted  by 
the  faculty  of  the  Department  of  Physics  and  Astronomy  at  the  University  of 
Mississippi  on  September  10,  1985.  The  work  reported  was  carried  out  at 
the  University  of  Mississippi's  Physical  Acoustics  Research  Laboratory  with 
support  from  the  Office  of  Naval  Research. 


Accession  For 

*NTIS 
DTIC  T ‘ 3 
Unannounced  □ 

Justification. 


Henry  E.  Bass 
Dissertation  Director 


LIST  OF  TABLES 


Page 

Table 


1.  Experimental  Rise  Times  and  Decay  Times  for  the  Positive 

Pressure  Amplitude  .  . . ...» . 35 

2.  Rise  Time  and  Decay  Time  for  the  Negative  Pressure 

Amplitude . •  . . 35 


iv 


LIST  OF  FIGURES 


Figure 

1 •  Block  Diagram  of  the  Experimental  Apparatus . 

2.  The  Microphone  Response . . . 

3.  Comparison  of  Eq.  (30)  and  Bass  and  Yan's  Result  .  .  .  . 

4.  A  Sequence  of  Images  of  Source  ......  . 

5.  Computed  Pressure  Evolution  at  High  Static  Pressure.  .  . 

6.  Computed  Waves  at  Intermediate  Pressures  . 

7.  Typical  Pressure  Waveform . 

8.  Change  in  Sign  of  Pressure  Amplitude  as  the  Gas  Pressure 

is  Decreased  . 

9.  The  Rise  Time  and  Decay  Time  ...  . 

10.  Partial  Vibrational  Energy  Level  Diagram  for  SF^  .  .  .  . 

11.  Computed  Curves  for  Sequence  of  Pressures . 

12.  Acoustic  Effects  at  Intermediate  Pressure  Curves  .  .  .  . 

13.  Fitting  of  Experimental  Data  with  Computed  Points.  .  .  . 

14.  Fitting  of  Experimental  Data  with  Computed  Curve  .  .  .  . 

15.  The  Effect  of  Spontaneous  Emission  at  Low  Pressure  .  .  . 


Page 

6 

10 

22 

24 

28 

30 

32 

33 

34 
37 

39 

40 

42 

43 

44 


v 


CHAPTER  1 


Introduction 

Alexander  Graham  Bell  (1880)^  observed  that  audible  sound  can  be 
produced  by  exposing  a  gas  in  a  constant  volume  container  to  intensity 
modulated  infra-red  radiation.  He  explained  that  part  of  the  radiation 
is  absorbed  by  the  gas  raising  the  temperature  and  hence  pressure  in  the 
constant-volume  enclosure.  A  periodic  modulation  of  the  incident 
intensity  thus  produces  an  alternating  pressure  component  which  becomes 
audible  for  proper  choices  of  the  experimental  parameters .  Bell  realized 
the  possibility  of  using  the  effect  to  investigate  the  absorption  spectra 
of  gases.  He  proposed  the  name  "spectrophone"  for  the  instrument 
constructed  for  this  purpose. 

2 

With  the  advent  of  molecular  theory  this  process  was  interpreted  in 
detail.  When  molecules  absorb  infra-red  radiation,  the  energy  taken  up 
appears  as  quanta  of  vibrational  energy.  A  molecule  may  lose  this  energy 
by  spontaneous  emission  of  radiation,  by  induced  emission,  or  by  transfer 
of  energy  to  thermal  energy  of  motion.  This  last  process  occurs  in 
collisions  when  the  vibrational  energy  is  transferred  to  translation.  An 
increase  in  translational  energy  means  the  gas  will  heat  and  reach  a  new 
temperature.  If  the  gas  is  illuminated  by  infra-red  radiation  which  is 
periodically  interrupted,  it  will  heat  and  cool,  and  a  periodic  change  of 
pressure  will  be  observed  in  a  constant  volume  system.  These  periodic 
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pressure  changes  lead  to  sound  emission  which  can  be  detected  by  means  of 
a  microphone. 

Interest  in  the  spectrophone  has  increased  after  the  experimental 

tools  applied  for  investigation  of  molecular  energy  transfer  in  gases 

shifted  from  pure  thermodynamic  ones  to  optical  methods.  Part  of  this 

interest  is  because  the  spectrophone  proves  to  be  a  powerful  tool  to 

measure  relaxation  rates.  Early  spectrophone  measurements  of  vibrational 

relaxation  used  a  light  source  modulated  to  be  nearly  sinusoidal  in 

amplitude  with  a  variable  frequency.  The  time  required  to  transfer 

energy  from  the  excited  vibrational  state  to  translation  was  measured  as 

either  frequency  dependent  phase  shift  or  signal  amplitude.  Filtered 

radiation  can  be  used  so  that  for  polyatomic  molecules  the  relaxation 

time  for  individual  vibrational  modes  can  be  studied  separately. 

3 

Gorelik  pointed  out  that  if  there  is  a  time  lag  between  the  absorption 
of  a  vibrational  quantum  and  its  degradation  to  thermal  energy,  it  should 
determine  a  phase  difference  between  the  modulation  of  the  incident 
radiation  and  the  emitted  sound  wave.  This  phase  shift  is  then  a 
function  of  the  relaxation  time  of  the  vibrational  level.  Also  the 
amplitude  of  the  sound  wave  is  a  function  of  the  relaxation  time  of  the 
absorbing  level. 

Modern  measurements  employ  either  modulated  or  pulsed  light. 

Measurements  of  vibrational  energy  transfer  using  an  amplitude  modulated 

4 

laser  beam  and  phase  sensitive  detection  have  been  most  common  and  have 
provided  a  great  deal  of  information  about  the  details  of  energy 
transfer.  The  pressure  response  following  illumination  depends  in 
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magnitude  (and  direction)  on  the  amount  of  energy  transferred  to  or  from 
translation.  This  feature  provides  the  spectrophone  with  a  unique 
capability  to  determine  energy  transfer  pathways. 

The  pulsed  spectrophone^  offers  certain  advantages  over  the 

continuously  modulated  configuration.  One  major  advantage  is  the  ability 

to  delay  interference  from  cell  resonances  by  varying  the  cell  length. 

An  additional  advantage  is  that  the  amplitude  and  time  evolution  are 

measured  as  separate  quantities  at  the  same  time  (during  a  single  shot) 

rather  than  measuring  the  phase  shift  between  modulated  excitation  and 

£ 

pressure  response  as  a  function  of  frequency.  This  feature  requires 
that  the  microphone  and  electronics  used  in  the  pulsed  spectrophone  have 
a  large  bandwidth. 

Interpretation  of  spectrophone  data  requires  accurate  relations 
between  the  microscopic  energy  transfer  processes  and  the  measured 
pressure  variation  as  well  as  a  description  of  other  macroscopic 
processes  which  affect  the  signal,  e.g.,  acoustic  propagation  and  thermal 
conduction.  Bauer^  has  developed  a  general  theory  for  relating  the 
microscopic  energy  transfer  processes  to  the  opto-acoustic  response  for  a 
gas  with  multiple  internal  energy  levels.  Basically,  using  his  approach, 
one  must  assume  a  set  of  energy  transfer  rates  and  a  reaction  scheme, 
then  compute  the  predicted  response.  The  resulting  pressure  evolution  is 
compared  to  experiment  and  the  rates  and  reaction  scheme  are  adjusted 
until  agreement  is  achieved.  There  is  no  assurance,  of  course,  that  the 
resultant  rates  and  reaction  scheme  are  unique.  This  is  true  of  other 
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methods  used  to  measure  vibrational  relaxation  times  as  well,  e.g.,  time 
resolved  laser  induced  fluorescence  and  ultrasonic  absorption  and  disper¬ 
sion.  One  major  strength  of  his  generalized  approach  is  that  the  same 
procedure  can  simultaneously  compute  sound  absorption  curves  and  laser 
fluorescence  decay  times  for  the  assumed  reaction  scheme.  By  comparing 
the  predictions  to  sound  absorption,  laser  fluorescence,  and  spectrophone 
data  at  the  same  time,  some  reaction  schemes  can  be  eliminated.  The 
reaction  scheme  which  gives  results  in  agreement  with  all  three  experi¬ 
ments  is  still  not,  necessarily,  unique,  but  at  least  it  will  be 

consistent  with  all  available  data. 

In  this  work  a  pulsed  spectrophone  has  been  used  to  determine  its 

applicability  to  measure  relaxation  rates  in  a  polyatomic  gas.  The  SF^ 

molecular  system  was  chosen  as  a  test  case  because  energy  transfer  in  SF^ 

has  been  thoroughly  studied  using  other  techniques  and  it  strongly 

2  8  9 

absorbs  CC^  laser  radiation.  Data  available  from  other  techniques  ’  ’ 
will  enable  us  to  eliminate  reaction  schemes  inconsistent  with  other 
techniques. 

The  objective  of  this  work  was  to  develop  a  method  to  solve  the 
kinetic  equations  for  the  spectrophone  response  including  spontaneous 
emission  and  acoustic  effects  as  a  function  of  time  and  to  test  the 
solution  experimentally.  This  objective  has  been  achieved  by  solving  the 


heat  equation  for  the  gas  inside  the  cell  by  numerical  integration. 


CHAPTER  2 


Description  of  Experiment 

A  block  diagram  of  the  experimental  apparatus  is  shown  in  Figure  1 . 
The  apparatus  consists  of  a  source  of  pulsed  1R  radiation,  a  test  cell,  a 
microphone  and  preamplifier,  signal  averaging,  and  gas  handling 
equipment. 

IR  Source 

The  IR  source  is  a  CO^  laser  consisting  of  a  2.5  cm  i.d.  glass  tube, 
125  cm  long  sealed  with  sodium  chloride  (index  of  refraction  1.46  at 
10. 6p)  windows  mounted  at  the  Brewster's  angle  of  56°,  and  two  aluminum 
coated  mirrors.  A  mixture  of  11%  CO^,  8%  and  81%  He  (Matheson,  Co. 
premixed)  flows  in  at  one  end.  An  electric  discharge  through  the  gas 
mixture  is  maintained  by  a  continuous  D.C.  power  supply  operating  at 
about  7kV  and  20  mA.  The  C02-N2_He  mixture  is  pumped  out  of  the  tube  at 
the  other  end.  The  continuous  flow  serves  to  cool  the  laser  gas  and 
helps  to  remove  impurities  and  vibrationally  excited  gas  which  might 
decrease  the  laser  efficiency.  Further  cooling  of  the  laser  gas  is 
provided  by  thermal  conduction  through  the  tube  walls  to  a  surrounding 
water  jacket.  Pressure  inside  the  laser  tube  was  generally  6-10  torr. 
The  CO2  laser  is  Q-switched  with  a  rotating  mirror  which  is  mounted  on  a 
D.C.  motor  to  deliver  pulses  approximately  1  ps  in  duration.  Mirror  M^ 
is  the  output  mirror  and  is  coated  for  20%  transmission  at  10.6  pm.  The 
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Block  Diagram  of  the  Experimental  Apparatus. 


laser  is  capable  of  delivering  approximately  1  W  of  10.6pm  radiation  when 
operated  in  a  continuous  mode. 

Associated  with  the  laser  is  a  triggering  circuit  used  to  initiate 
the  electronics.  This  circuit  consists  of  a  1  mW  He/Ne  laser  and  a  box 
containing  a  photodiode  cell,  amplifier  chip,  and  a  logic  circuit.  The 
power  to  the  circuit  was  provided  by  a  9V  battery.  The  He/Ne  laser  is 
directed  into  Mirror  Mj  so  that  just  before  rotates  into  a  position  to 
close  the  laser  cavity,  the  He/Ne  beam  strikes  the  photodiode  which 
provides  a  pulse.  After  amplifying  the  pulse  with  an  amplifier  chip  and 
shaping  it  by  the  logic  circuit  to  give  a  5V  square  wave,  it  was  used  to 
trigger  the  signal  averager. 

Test  Cell 

The  body  of  the  cell  is  a  tube  mounted  so  that  its  axis  is  coincident 
with  the  laser  tube  axis,  and  hence  the  laser  beam  axis.  The  character¬ 
istics  of  the  test  cell  are  a  strong  function  of  the  particular  system  to 
be  studied.  Cylindrical  geometry  seems  natural  due  to  the  approximately 
cylindrical  beam  shape.  The  length  of  the  test  cell  is  chosen  to  be  50 
cm  long.  Selection  of  the  cell  radius  is  influenced  by  several  interde¬ 
pendent  considerations.  Across  the  cell  radius,  the  gas  is  illuminated 
uniformally  and  the  pressure  evolution  occurs  simultaneously.  The 
microphone,  however,  is  located  in  the  side  of  the  test  cell  so  it 
responds  to  pressure  changes  at  the  cell  center  after  a  time  where  r  is 
the  cell  radius  and  c  is  the  speed  of  sound  in  the  gas. 
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The  resulting  pressure  response,  then,  will  be  averaged  over  a  time 
^  and  any  processes  with  a  shorter  time  constant  will  not  be  distin¬ 
guished  clearly.  For  SFg,  c  s  100  m/sec  so  in  order  to  observe  a  20  psec 
relaxation  time,  r  must  be  less  than  2  mm.  Transit  time  considerations 
argue  for  a  small  r.  When  the  radius  is  decreased,  however,  thermal 

conduction  becomes  more  important.  The  characteristic  time  for  thermal 

2 

conduction  x  _  ,  decreases  as  r  .  If  x  _  is  much  less  than  the  slowest 
tc’  tc 

relaxation  process  of  interest,  that  slow  process  will  be  difficult  to 
observe.  The  final  consideration  is  the  magnitude  of  the  pressure 
response.  In  order  to  avoid  multiphoton  processes,  the  laser  intensity 
cannot  be  increased  indefinitely.  As  the  cell  cross  section  decreases, 
then,  the  signal  amplitude  will  decrease. 

We  selected  a  50  cm  long  test  cell  with  radius  of  2  mm.  The  cell  is 
approximately  cylindrical  with  salt  flats  on  both  ends.  The  microphone 
is  mounted  in  the  side  of  the  test  cell  near  the  entrance  window.  The 
microphone  extends  slightly  into  the  2  mm  cylinder  so  the  effective 
radius  at  the  microphone  location  is  less  than  2  mm,  perhaps  1  mm  from 
beam  center  to  the  nearest  part  of  the  microphone  face.  Gas  inlets  are 
provided  near  each  end  of  the  cell  so  that  the  cell  can  be  flushed  with 
the  gas  of  interest.  Also  connected  directly  to  the  test  cell  are  the 
pressure  transducers. 

Microphone 

The  microphone  is  very  important  foe  the  measurements.  The  response 
of  the  microphone  can  affect  the  recorded  pressure  response.  The 
microphone  response,  at  low  frequencies,  is  determined  by  leakage  of  air 
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from  front  to  back  of  the  diaphragm.  At  high  frequencies,  the  diaphragm 
begins  to  oscillate  in  a  resonant  mode.  For  pulsed  measurements,  the 
microphone  must  be  broad-band  with  an  upper  frequency  limit  fu»l/(2uT) 
where  t  is  the  fastest  relaxation  process  of  interest.  For  our  purpose 
f  should  be  greater  than  20  kHz.  This  condition  is  easily  achieved  with 
modern  capacitance  microphones. 

A  standard  commercial  capacitance  microphone  has  an  air  gap  between  a 
backing  plate  and  a  metal  diaphragm  which  forms  the  pressure  sensitive 
capacitance.  These  microphones  rely  on  this  air  gap  to  damp  the 
diaphragm  motion  and,  as  a  result,  quite  often  ring  when  hit  by  a  short 
pulse  at  low  pressure.  For  this  reason,  we  used  a  solid  dielectric 
capacitance  microphone  with  an  upper  frequency  limit  of  approximately  100 
kHz.** 

The  low  frequency  microphone  response  was  measured  using  a  substitu¬ 
tion  technique  with  a  B&K  4149  microphone  as  the  standard.  The  speaker 
used  as  a  source  consisted  of  a  16  inch  diameter  driver  mounted  in  a  725 
liter  box  giving  a  frequency  response  of  50-5000  Hz.  The  electrical 
source  was  taken  from  an  HP  356 1A  dynamic  signal  analyzer.  The  analyzer 
generates  band  limited  noise  then  Fourier  analyzes  the  input  signal.  The 
signal  from  the  reference  microphone  is  stored  and  subtracted  from  the 
signal  from  the  microphone  under  test.  Results  of  the  calibration  proce¬ 
dure  are  shown  as  Figure  2.  Note  that  no  effort  was  made  to  define  the 
high  frequency  cut  off  of  the  test  microphone.  We  did  determine  that 
this  cut  off  was  above  that  of  our  reference  microphone  (30  kHz).  Since 
the  microphone  relies  on  gas  in  pockets  between  the  dielectric  membrane 


Figure  2.  The  Microphone  Response 
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and  the  backing  plate,  the  sensitivity  decreases  slightly  with  lower 
pressure,  but  the  loss  in  sensitivity  is  much  more  acceptable  than 
ringing. 

The  microphone  signal  was  preamplified  first  with  a  General  Radio 

type  1560-42  preamplifier  specifically  designed  to  amplify  the  output 

from  a  capacitive  source.  The  preamplifier  offers  a  gain  of  either  XI  or 

X10.  A  self-contained  oscillator/rectifier  circuit  supplies  the 

polarizing  voltage  (+200V)  for  air-condenser  microphones.  The  signal  was 

preamplified  again  with  an  EG&G  PARC  Model  113  preamplifier.  Amplifiers 

were  flat  within  3  dB  from  below  7  Hz  to  above  200  kHz.  Total  gain  was 
4 

typically  10  . 

Signal  Averager 

The  preamplified  microphone  signal  was  averaged  with  a 

Honeywell/ Saicor  averager.  Its  function  when  operated  in  the  signal 

enhancement  mode  is  to  provide  detection  and  enhancement  of  a  periodic  or 

periodically  stimulated  signal  buried  in  noise  or  other  disturbance 

through  an  averaging  procedure  which  preserves  the  wave  shape  of  the 

signal  while  averaging  out  non-coherent  noise.  The  waveform  is  divided 

into  400  points  with  resolution  (or  spacing  between  points)  selectable 

9  17 

from  0.2  psec  to  1  sec.  From  2  to  2  (in  binary  steps)  successive 
cycles  can  be  linearly  averaged  in  addition  to  the  continuous  operation 
mode.  The  averager  was  triggered  externally  by  the  He/Ne  laser  reflected 
from  the  rotating  mirror  which  Q-switched  the  C0^  laser.  The  output  S/N 
is  typically  less  than  one  so  the  signal  cannot  be  viewed  on  an  oscillo¬ 
scope  except  at  high  pressure.  The  signal  averager  was  operated  either 
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In  the  clipped  or  full  mode  depending  upon  S/N  ratio.  The  output  from 
the  averager  was  sent  to  an  X-Y  recorder  for  a  permanent  record. 

Gas  Handling  Equipment 

The  SF^  gas  used  was  supplied  by  Matheson  with  a  stated  purity  of 
99.982.  To  keep  water  vapor  from  leakage  and  outgasslng  low,  the  gas 
cell  was  evacuated  for  several  days  prior  to  measurements.  The  combined 
leak/outgassing  rate  was  0.3  mtorr/min.  Pressures  were  measured  with  MKS 
Baratron  tranducers  type  220.  The  transducers  were  factory  calibrated 
and  certified  for  0.1  percent  accuracy. 

The  experiment  consists  of  turning  on  the  laser  and  the  rotating 
mirror  to  Q-switch  the  laser  and  adjusting  the  mirrors  to  get  radiation 
of  maximum  energy.  Energy  was  measured  with  a  pyroelectric  energy  meter 
to  be  “20  uJ/pulse.  Then  the  cell  was  filled  with  SF^  gas  to  the  maximum 
desired  pressure  of  measurement  “1000  mtorr.  Next  the  cell  was  illumi¬ 
nated  with  TR  radiation  from  the  laser.  The  test  cell  had  IR  windows  on 
both  ends  to  allow  the  laser  beam  to  enter  and  exit. 

The  acoustic  signal  was  picked  up  by  the  microphone  which  was 
mounted  in  the  cell  wall  near  the  entrance.  At  the  same  time  the  He/Ne 
laser  was  turned  on  to  initiate  the  triggering  of  the  signal  averager. 
After  the  signal  had  passed  through  the  preamplifier  and  signal  averager, 
it  was  recorded  on  the  X-Y  recorder.  The  pressure  in  the  test  cell  was 
then  lowered  to  the  next  desired  pressure  and  another  run  was  made.  To 
avoid  accumulation  of  water  vapor,  the  test  cell  was  evacuated  every 
three  runs.  This  continued  until  the  desired  pressure  range  had  been 


covered 


CHAPTER  3 
Theory 

In  1972,  Bauer^  published  a  theory  for  the  spectrophone  response 
results  from  collisional  energy  transfer  using  a  matrix  notation  which  Is 
well  adapted  to  numerical  calculation.  Thermodynamic  conditions  are 
assumed  to  be  adiabatic.  Bauer  assumed  the  system  was  perturbed  by 
optical  transitions  then  the  resulting  collisional  transitions  were 
considered.  The  collisional  change  of  the  population  numbers  determines 
the  amount  of  energy  exchanged  between  the  translational  and  internal 
degrees  of  freedom.  When  applying  the  theory  to  a  multilevel  gas,  Bauer 
pointed  out  that  application  of  linear  algebra  allows  one  to  represent 
the  matrices  by  their  eigen-properties  which  leads  to  a  decomposition  of 
the  optoacoustic  response  into  single  relaxation  terms. 

The  main  steps  of  Bauer's  development  will  be  used  here  to  derive  the 
temperature  response  resulting  from  absorbed  radiation  (for  more  infor¬ 
mation  cf.  Ref.  7).  Bauer  considered  multiple  optical  transitions  as 
well  as  collisional  transitions. 

Collision  processes  refer  to  transferring  vibrational  energy  from  one 
mode  to  another  or  to  translation.  Thermal  conduction  is  also  a  result 
of  energy  transfer  during  collisions  but  for  the  purpose  of  this  disser¬ 
tation  thermal  conduction  will  not  be  included  in  the  term  collisional 
processes.  Thermal  conduction  will  be  treated  separately. 

Optical  Transition 

Consider  one  optical  transition  which  takes  a  molecule  of  type  M  in 
state  k  (denoted  M^)  to  state  i  by  addition  of  a  quantum  of  radiation 
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with  frequency  The  reaction  can  be  written  as 

Mi  *  Vhvik  (1) 

The  rate  of  change  of  occupation  numbers  (n^  for  population  of  and  ^ 
for  population  of  M^)  due  to  absorption  and  stimulated  emission  can  be 
written  as 

"i  "  vik(BkiVBikni)  (2) 

where  B^  and  are  the  Einstein  coefficients  for  absorption  and  stimu¬ 
lated  emission  respectively.  In  Eq.  (2),  spontaneous  emission  has  been 
neglected.  With  this  assumption,  the  rate  of  change  of  occupation  number 
due  to  optical  transitions  can  be  written  using  matrix  notation  as 

i.Pt  ■  v  <3) 


Colllsional  Transitions 

Bauer  introduced  inelastic  binary  collision  processes  of  the  type 

k' 

a 

Mi+Mj  z  VMi  (*> 

where  . may  be  any  energy  level  of  the  gas  or  mixture  considered 

and  k&  and  k^  are  the  forward  and  backward  rates  in  units  of  sec  *atra 
The  net  rate  of  the  reaction  is  given  by  the  difference  between  forward 
and  backward  rate. 

The  rate  of  the  a1*1  reaction  is 

«« *  np(Vixj^xkxt)  (5) 

where  £  is  the  progress  variable  for  reaction  a,  p  is  the  pressure  in 
a 

atmospheres,  n  is  the  total  number  of  molecules  and  x^.  .  .x^  are  the 

mole  fractions  (x^  ■  n^/n). 


C  Is  then  expanded  about  a  steady-state  value  assuming  a  small 


perturbation.  Using  the  mass  action  law, 

31n(pk  )  31n(pk') 

1  “  npk°x°x°  {[ - g— - W - ](T-T°) 

3 In (x  x  )  31n(x  x  ) 

+  I[ — - -n°)>  (6) 

m  3n  an  m  m 

m  m 

where  superscripts  o  refer  to  the  value  at  the  expansion  point.  Bauer 

related  the  first  bracket  to  the  reaction  heat  AH  and  used  Van  t'Hoff's 

a 

equation  to  write 


—  -  AH  /RT2. 


The  second  bracket  can  be  rewritten  as 


31n(x  x  )  3ln(x  x  ) 

rr  i  i  k**i/  0\ 

£[ - — 4 - sn - 1  (n  -n  ) 

m  3n  3n  m  ra 

m  m 

o  o 


-1  r  ni-ni  .  n1-ni 

n  [ - +  — 1 1 — ■*- - ] . 

Xi  Xi  \  XZ 


For  multilevel  systems  a  matrix  notation  is  introduced.  Let  be  the 
vector  of  ?a;  ri  be  the  vector  of  n^-n°;  ]i  be  c^e  diagonal  matrix  (k^x^x^ ) 
and  3t  be  the  diagonal  matrix  of  all  mole  fractions. 

Further,  Bauer  constructed  a  matrix  _v  with  columns  corresponding  to 
an  individual  level,  such  that  (v)  counts  the  number  of  times  the 
energy  level  i  occurs  more  often  on  the  right  hand  side  of  the  reaction 
equation  a  than  on  the  left  hand  side.  This  matrix  allows  one  to  express 


all  AH  in  terms  of  the  energies  of  one  mole  of  each  level,  accumulated 
a 

in  the  vector  u,  as  \3u.  Equations  6-8  now  read  as 
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^  =  -Pkv  ljfln-(^)u(T-T°)]. 


(9) 


Multiplying  £  by  \>,  the  rate  of  change  of  the  mole  numbers  due  to  colli¬ 
sion  processes  is  found  to  be 


^coll  "  "P^k5  u(T-T°)] 


(10) 


The  perturbation  term  n^^  generates  a  temperature  variation  with 


the  rate 


T  -  -utl  /(nC°°)  (11) 

— coil  v 

•  00 

where  -un  is  the  rate  of  energy  flow  into  translation,  and  Cv  is  the 

molar  heat  capacity  of  translation  and  all  degrees  of  freedom  which  share 
their  energy  instantly  with  translation. 


Spontaneous  Emission 

Bauer  neglected  the  effects  of  spontaneous  emission.  Spontaneous 
emission  has  two  experimentally  observable  effects.  At  high  pressure, 
since  the  effect  of  thermal  conduction  decreases  with  increasing 
pressure,  spontaneous  emission  is  a  mechanism  for  transferring  energy 
supplied  to  the  gas  by  the  laser  to  the  cell  walls.  This  effect  will 
parallel  thermal  conduction  and  will  be  important  for  high  gas  pressures 
and  long  times.  At  low  pressure,  spontaneous  emission  will  offer  a 
parallel  relaxation  path  to  the  slow  collisional  processes.  As  a  result, 
the  observed  relaxation  times  will  be  a  combination  of  collisional, 
spontaneous  emission,  and  wall  deexcitation. 
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The  rate  of  change  of  occupation  number  due  to  spontaneous  emission 


is 


n  *  -An 
s  — 


(12) 


where  A  is  the  inverse  of  the  radiative  lifetime.  Including  spontaneous 
emission,  Equation  (3)  becomes 


n  =*  n ,+n  =  -An+n, . 

—opt  — i  — s -1 


(13) 


The  total  rate  of  change  of  occupation  number  is  now 


n  -  n  +  n 
—  —opt  —coll 


(14) 


Inserting  Eq.  (13)  into  Eq.  (14)  yields: 


*  •  •  . 
n-n .  =  n  , , -An . 
- i  —coll  — 


(15) 


Eq.  (15)  can  be  solved  numerically.  At  time  t  *  0,  the  change  in  temper- 

o  •  • 

ature  T-T  *  0,  then  n,  and  n  ,,  can  be  determined  by  using  Eqs.  (2)  and 

— I  —coll 

(10)  respectively,  and  then  n  can  be  computed  from  Eq.  (15).  At  a  later 

+0 

“ “•  • 

time,  t  =  At,  a  new  value  for  n  =  f  n  dt+nAt  is  computed.  Also  a  new 
o  "  -0-Pt  " 

value  for  T-T  is  determined  by  using  Eq.  (11).  These  new  values  for  T- 


T  and  n_  are  used  in  Eq.  (10)  to  compute  a  new  value  for  n^  .  Eq*  (15) 
is  then  used  to  determine  a  new  value  for  n. 

By  repeating  this  procedure  again  and  again  for  small  time  changes, 
the  entire  temperature  evolution  can  be  computed.  The  pressure  change  is 
then  computed  from  the  temperature  change  using  the  ideal  gas  law. 
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Thermal  Conduction 

In  solving  the  kinetic  equations  for  the  spectrophone  response,  Bauer 
assumed  that  thermal  conduction  to  the  cell  walls  is  very  slow  compared 
to  the  rate  of  collisional  deactivation.  This  Is  true  at  high 
pressures.  At  low  pressures  of  interest  here,  however,  thermal  conduc¬ 
tion  can  dominate  the  measured  response. 

The  equation  for  heat  flow  in  the  test  cell  when  there  are  sources  of 

10 

heat  within  it  can  be  written  as 


-  *<-■£&  *  4  .  9(r>t) 


where  9  is  the  difference  between  the  gas  and  tube  temperatures  (the  tube 
temperature  is  assumed  to  be  constant).  The  rate  of  heat  generation, 
H(r,t),  is  assumed  to  be  a  product  of  space  and  time  dependent  parts, 


H(r,t)  »  f (r)g(t)  so  0(r,t)  =  -^-H(r,t)  =  ^-f(r)g(t)  (17) 

where  k  is  the  "thermal  diffusivity"  and  k  is  the  thermal  conductivity 


(k  »  kC^p).  The  function  g(t)  is  the  heat  generated  by  relaxation 


processes  which  is  equal  to 


3u(t) 


Prior  to  proceeding  further  with  the  solution  of  Equation  (16),  a 
form  for  f(r)  must  be  assumed.  We  will  assume  a  uniform  radial  excita¬ 
tion  though  some  type  of  a  Gaussian  distribution  across  the  cell  radius 
might  be  more  appropriate. 

For  uniform  excitation  f(r)  «  1  and 
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Equation  (16)  can  then  be  written  as 


£[ecr,t)  -£.<«]  -  11211^1. 


(19) 


Equation  (19)  can  be  solved  by  separation  of  variables.  Let 


0(r,t)  =  R(r)Q(t). 


(20) 


Substituting  0(r,t)  from  Equation  (20)  into  Equation  (19)  and  dividing  by 
R(r)Q(t)  yields 


cR(r)Q (t)  3ttR(r)Q(t)  “  ku(t)]  R(r)Q(t)  ^3r5'R(r)C!^t) 


)Q(t 

+  ^R{r)Q(t)]  . 


(21) 


By  equating  Equation  (21)  to  -a  ,  the  radial  equation  (the  right  hand 

side),  has  a  solution,  R  =  AJ  (or),  where  J  is  Bessel's  function  and  A 

o  o 

is  a  normalization  constant  to  be  determined  from  the  boundary 
conditions. 

The  boundary  conditions  require  that  9=0  when  r  =  the  cell  radius 


a,  hence  J  (a  a)  =■  0  and  a,  ,  a 
o  in  l  L 

function.  The  radial  solution  is,  then, 


a  are  zeros  of  the  Bessel 
m 


R  (O 

m 


l  A  J  (a  r). 
m  m  o  to 


(22) 


A  can  be  determined^  to  be  21 


m  aa  J, (aa) 
tn  1 


so 


J  (a  r) 

R  (r)  *  21  —  -7 — r-. 
ra  aa  -.(aa) 

m  i 


(23) 
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Substituting  for  R  into  Eq.  (20)  yields 
m 

J0(a  r) 

0(r»t)  “  2&  i3r°j-Cct-a)Q(t)- 

m  1  m 


Taking  the  spatial  average  of  Eq.  (24)  by  integrating  across  the  cell 
radius , 


0(t)  ■  pr  f  ®(r ,t)dr  =  4Q(t)£  £?£? 

0  m  m 


where 


Q(t)  --58(c)  [ 


The  left  hand  side  of  Eq.  (21)  can  be  written  as 


ft  lQ(t)£  \(t)  -  !“Wi  - 


Averaging  over  r. 


[^t4t-Q(c)  -  -p(t)]  =  -  |r  Q(t) 
m 


Substituting  for  Q(t)  into  Eq.  (28) 


k  du  ,  .  K0  ,  N 

-kS(t)--p(c) 


If  we  substitute  for  k  and  — ,  Equation  (29)  becomes 


M  i  Pvk\3  fx""^n  -  (■  n  >  fT  T°>  1  ^RT 

dt  nC”  — -  T  ^  Pr“  2-1 


PC"  a2  ma2a  2 
v  m 


Eq.  (30)  consists  of  two  terras.  The  first  one  is  due  to  collisional 


processes  given  earlier  in  Eq.  (11)  and  it  is  directly  proportional  to 


the  pressure.  The  second  one  is  due  to  thermal  conduction  and  it  i9 
inversely  proportional  to  the  pressure. 

It  is  clear  from  Eq.  (30)  that  at  high  pressure  the  thermal  conduc¬ 
tion  term  can  be  neglected  and  the  rate  of  change  of  temperature  is 
mainly  due  to  collisional  processes,  while  at  low  pressure  thermal 
conduction  can  not  be  neglected.  Eq.  (30)  can  be  solved  numerically. 

Bass  and  Yan^  have  modified  Bauer's  theory  by  including  thermal 
conduction  and  wall  collisions  in  the  theory.  They  also  used  Bauer's 
matrix  notation  to  solve  for  the  kinetic  equations  of  the  spectrophone. 
Our  solution  of  Eq.  (30)  using  numerical  integration  is  compared  to  the 
Bass  and  Yan  result  in  Figure  3 . 

Effect  of  Wall  Collisions 

12 

According  to  Margottin-Machou,  Dogonette,  and  Henry  ,  the  effective 

relaxation  rate  &(*  — — )  can  be  written  as 

Teff 

g  -  k  +k  +iv  (31) 

c  w  p 

where  k  is  the  collisional  deactivation  rate,  k  is  the  rate  of  deexci- 
c  w 

tation  by  the  walls  and  is  the  rate  of  deactivation  by  spontaneous 
emission  (including  trapping).  They  further  show  that  for  a  long  narrow 
tube, 

k  -  p2(D/R2)  (32) 

w 

where  D  is  the  diffusion  coefficient,  R  is  the  cell  radius,  and  p  is  a 

solution  of  pJ.(p)  *  hJ  (p),  where  h  ■  (Rv/2D)(a/2-a) ,  v  is  the 

l  o 


Comparison  of  Eq.  (30)  (Lower  Curves)  and  Bass  and  Yan's  Result  (Upper  Cu 
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mean  speed  of  molecules,  and  a  is  the  probability  of  deexcitation  during 

a  collision  with  the  cell  walls.  The  ratio  R/h  is  the  minimum  distance  a 

molecule  must  travel  close  to  a  wall  before  being  deexcited  sA/a  where  A 

is  the  mean  free  path.  The  diffusion  coefficent  D  is  pressure  dependent 

so  we  will  write  D  **  DP. 

o 

Margottin-Machou,  Dogonette,  and  Henry  found  that  for  pressures  of  1 
Torr  and  above,  p  approaches  a  limiting  value  of  2.4.  For  pressures  of 
interest  in  this  work  p  varies  from  2.4  -  3.6.  Changing  from  2.4  to  3.6 
has  no  significant  effect  on  the  shape  of  the  pressure  wave,  so  Eq.  (31) 
was  written  as 

k  =*  (2.4)2D  /R2P  (33) 

w  o 

where  D  is  the  kinetic  self-diffusion  coefficient, 
o 

Acoustic  Effect 

As  the  incoming  radiation  propagates  down  the  axis  of  the  tube, 

absorption  of  radiation  by  the  gas  decreases  the  intensity  in  accordance 

with  Beer's  Law,  I  =  I  e  .  Since  the  absorbed  energy  in  any  small 

o 

section  of  the  tube  differs  from  that  in  the  next  section,  the  resulting 
temperature  and  pressure  rise  will  decrease  down  the  tube  axis.  A  pres¬ 
sure  gradient  down  the  tube  will,  of  course,  give  rise  to  a  net  force 
causing  gas  near  the  entrance  window  to  expand  into  the  lower  pressure 
regions  further  from  the  entrance  window.  This  expansion  is  observed  as 
a  longitudinally  propagating  acoustic  wave  (superimposed  upon  a  uniform 
pressure  change).  For  the  acoustic  wave,  the  effect  of  the  boundaries  is 
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important.  Boundaries  cause  reflections  which  must  be  added  to  the  wave 
developed  by  the  source  to  give  the  total  pressure.  To  determine 
the  total  pressure  at  any  point,  we  will  use  a  Green's  function  approach 
and  the  method  of  images. 

13 

According  to  Morse  and  Feshbach  the  solution  to  the  one-dimensional 
scalar  wave  equation  is 
t  x 

p(x,t)  3  ^0dto  /gdx^x.t  |xo,to)T(xo,tQ)  (34) 

where  xq  represents  the  source  point,  tQ  represents  the  time  when  the 

source  is  applied,  x  represents  the  observation  point,  where  the  field  is 

measured,  t  represents  a  time  later  than  when  the  source  was  applied  at 

x  ,  p  is  the  presssure  at  the  observation  point,  T  is  the  temperature 
o 

which  describes  the  source  density  as  a  function  of  space  and  time,  and  G 
is  the  Green's  function  describing  the  effect  of  the  source  as  it 
propagates  away  from  xq  in  the  course  of  time. 

To  apply  the  method  of  images  consider  the  test  cell  as  a  cylindrical 
tube  of  length  L  as  shown  in  the  following  diagram. 


Microphone 

x  =  0  |  +  \ _ x  =  L 


1 

» 

I 


-X 


x  =  2L 

~2l 

•  I 

_  J 


Figure  4.  A  Sequence  of  Images  of  Source. 
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The  method  of  images  is  applied  in  a  stepwise  fashion,  by  considering 

each  boundary  x  =  0  and  x  *  L  in  turn.  Disregarding  the  x  ■  L  boundary, 

an  image  is  introduced  at  x  =  -xq.  Disregarding  x  *  0,  an  image  is 

introduced  at  x  »  2L  -  x  .  This  process  could  continue  indefinitely, 

o 

leading  to  an  infinite  number  of  images,  since  there  are  an  infinite 
number  of  reflections  for  any  ray  started  at  the  source.  Kach  image 
corresponds  to  one  of  these  reflections. 

The  number  of  images  required  will  depend  on  the  distance  travelled 
by  sound  inside  the  gas  sample  during  the  time  of  interest.  For  SF^,  the 
sound  travels  40  cm  in  the  4  msec  typically  of  interest  in  our  experi¬ 
ments.  Since  the  length  of  the  test  cell  is  50  cm,  this  Implies  that 
only  one  image  at  x  =  -xq  is  needed. 

13 

The  one  dimensional  Greeu.'s  function  for  our  problem  is  given  by 

I  x-x  I 

G(x, 1 1 x  , t  )  =  2Ctt [  1-u( - -  -  (t-t  ))]  (35) 

o  o  u  o 
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After  a  time  t  ■  up  to  t  ■  (L-x)/C,  there  is  one  image  to  consider 

located  at  -x  so 
o 

t- ~  x  =x+C (t-t  )  t  x  =x+C (t-t  ) 

f  C  (•  o  o  f  to  o 

p(x,t)  «  2tt C  dt  dx  T(x  ,t  )  +  2irC  dt  dx  T(x  ,t  ) 

JoJ  oo’  o  J°J  ooo 

0  Xo=x-C(t-tQ)  t=|  xq=0 

t  0 

+  2ttC  dt  dx  T(x  ,t  ).  (38) 

Jo  J  o  o  o 

t“f  Xo“X_C(t_to) 

To  solve  for  p(x,t)  we  need  to  determine  T(x  ,t  ).  T(x,t)  = 

*  '  O'  o  0^0 

-px0 

e  T  (t  )  where  p  is  the  absorption  coefficient  for  the  gas.  In  order 
o 

to  determine  T(t  )  the  standard  attenuation  mechanisms  must  be  considered, 
o 

These  losses  are  caused  by  viscosity,  thermal  conductivity,  mass 
diffusion,  and  vibrational  and  rotational  relaxation.  In  addition,  there 
are  losses  associated  with  sound  propagation  inside  the  tube.  When  sound 
propagates  through  polyatomic  gases  the  compressions  and  rarefactions 
involve  a  redistribution  of  energy  between  the  various  modes  of  excita¬ 
tion;  this  includes  vibrational  and  rotational  modes.  At  sufficiently 
high  frequencies,  the  period  of  fluctuation  becomes  shorter  than  the  time 
required  for  the  redistribution  of  energy  into  various  modes.  When  this 
occurs  the  compressions  and  rarefactions  are  no  longer  reversible  and 

attenuation  results.  Since  rotational  energy  levels  are  closely  spaced 

14 

near  room  temperature  for  SF^,  a  single  relaxation  time  can  be  used  to 

describe  rotational  relaxation.  For  vibrational  attenuation,  the 

15  16 

coupling  between  v-t  and  v-v  processes  complicates  the  problem  ’ 
greatly;  however,  SF,  has  properties  which  simplify  the  problem.  SF,  has 
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many  vibrational  levels  which  are  excited  even  at  low  temperatures  ~300K 
but,  the  lowest  level  deexcites  much  slower  than  the  others.  This  lowest 
level  gating  makes  it  possible  to  characterize  SF,  vibrational  relaxation 

D 

with  one  relaxation  tirae^. 

The  attenuation  coefficient  was  determined  by  using  a  computer 
program  developed  by  Shields  which  is  described  in  Reference  17. 

T(t^)  can  be  determined  by  the  following  steps: 

1.  Solve  Eq.  (30)  numerically,  T(t)  »  0(t)  can  be  determined. 

2.  Fourier  transform  T(t)  to  get  T(w). 

3.  Compute  the  attenuation  coefficient  a  as  indicated  earlier  for 
all  frequencies  of  interest. 

4.  Compute  T  (w)  *  T(w)10 

co 

5.  Inverse  T  (w)  to  get  T(t  ),  and  then  substitute  the  value  of 

co  o 

T(tQ)  into  Eq.  (38)  and  solve  for  p(x,t). 

Test  of  the  Computer  Program 

To  solve  Eq.  30  numerically,  a  computer  program  was  developed.  As  a 
first  test,  the  output  of  the  program  was  examined  for  pressures  of  1 
torr  and  2  torr.  At  these  pressures,  collisional  processes  dominate  the 
results.  Since  the  rate  of  collisional  deexictation  (one  over  the 
relaxation  time)  is  directly  proportional  to  pressure,  the  two  curves 
should  overlap  when  the  2  torr  curve  is  plotted  on  a  scale  covering  1/2 
the  time.  The  result  of  this  test  is  shown  in  Figure  5.  There  is  a 
slight  difference  due  to  the  finite  time  increment  used  for  the  numerical 
Integration.  This  difference  is  not  considered  significant. 


Figure  5.  Computed  Waves  at  High  Pressures.  —  for  pressure  =  1  torr;  —  for  pressure  =  2  torr 
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Figure  6  shows  computed  curves  at  intermediate  pressures  (0.8  and  0.4 
torr)  where  thermal  conduction  is  important  as  well  as  collisional  energy 
transfer.  The  long  term  return  of  the  pressure  back  to  equilibrium 
reflects  thermal  conduction  to  the  tube  walls.  From  Fig.  6,  the  thermal 
conduction  decay  time  for  a  pressure  of  0.8  torr  is  two  times  larger  than 
that  for  a  pressure  of  0.4  torr.  This  agrees  with  the  theory  which 
suggests  that  thermal  conduction  decay  time  (one  over  the  rate)  is 
linearly  proportional  to  the  pressure. ^ 


Acoustic  Pressure  (Arbitrary  Units)  Acoustic  Pressure  (Arbitrary  Units) 
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Time  (msec) 


Figure  6.  Computed  Waves  at  Intermediate  Pressures. 
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CHAPTER  4 


Experimental  Result 

A  typical  pressure  waveform  obtained  in  this  experiment  is  shown  in 
Figure  7.  This  figure  shows  the  variation  of  pressure  amplitude  (in 
arbitrary  units)  versus  time.  Energy  deposited  in  the  gas  by  a  short 
laser  burst  appears  first  as  a  small  perturbation  in  the  population  of 
vibrationally  excited  states.  Following  a  sufficient  number  of 
collisions,  the  excess  energy  will  redistribute  over  internal  and 
translational  modes.  Changes  in  temperature  and  pressure  will  result  as 
the  energy  cascades  down  through  the  manifold  of  vibrational  states.  The 
resulting  pressure  evolution  will  appear  as  shown  in  Figure  7.  Figure 
7  also  shows  that  the  pressure  decreases  (positive  pressure  is  down) 
for  very  short  times.  Since  there  is  a  net  gain  of  energy  by  the  gas,  on 
a  longer  time  scale,  the  pressure  will  go  positive.  For  even  longer 
times,  the  pressure  will  depend  on  acoustic  effects,  spontaneous 
emission,  and  thermal  conduction  to  the  cell  walls  which  take  the  system 
back  to  equilibrium. 

As  the  gas  pressure  is  decreased  the  amplitude  of  the  pressure  wave 
becomes  smaller  down  to  200  mTorr.  At  this  point,  the  sign  of  the 
pressure  change  reverses  (becomes  negative)  and  the  amplitude  begins  to 
Increase  as  the  pressure  is  decreased  further  as  shown  in  Fig.  8. 
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Sign  of  Pressure  Amplitude  as  the  Gas  Pressure  is  Decreas 
down 
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The  form  of  the  differential  equations  in  the  preceeding  Chapter 
suggests  that  the  pressure  will  vary  as  a  sum  of  exponential  increases 
and  decays,  i.e., 

p  a  p  e-t^T 

oe 

where  Pq  is  the  equilibrium  pressure  and  t  is  some  effective  decay  time 
which  will  change  with  time.  By  plotting  the  pressure  amplitude  on  a 
semilog  scale  versus  time,  rise  times  and  decay  times  for  specific  time 
intervals  can  be  determined  at  different  values  of  pressures  as  shown  in 
Fig.  9.  Table  I  shows  the  experimental  rise  times  and  decay  times  for 
the  positive  pressure  amplitude  (intermediate  pressure  regime).  It  can 
be  seen  from  this  table  that  the  relaxation  times  are  only  weakly 
dependent  on  pressure.  Table  2  shows  the  rise  time  and  decay  time  for 
the  negative  pressure  amplitude  (below  200  mTorr).  At  these  pressures 
both  times  are  pressure  dependent. 


Figure  9.  The  Rise  Time  and  Decay  Time  (Positive  Pressure  is  Down) 
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TABLE  I 


Pressure  (mTorr) 

Risetime  (msec) 

Decay  time  (msec) 

1000 

0.078 

.295 

800 

0.078 

.280 

600 

0.078 

.280 

400 

0.080 

.270 

300 

0.082 

.150 

TABLE  II 

• 

Pressure  (mTorr) 

Risetime  (msec) 

Decay  time  (msec) 

200 

.280 

1.98 

160 

.132 

1.84 

100 

.112 

1.44 

70 

.096 

1.37 

30 


094 


1.18 


CHAPTER  5 

Interpretation  and  Conclusions 

The  SF^  molecule  has  15  vibrational  degrees  of  freedom  and  six 

18  19 

fundamental  modes.  It  strongly  absorbs  10.6  pm  CO2  laser  radiation 

because  most  of  the  absorbing  molecules  are  in  the  ground  or  lowest 

excited  vibrational  state.  To  interpret  the  experimental  data  of  Figs.  7 

and  8  first  consider  the  energy  level  diagram  of  SF^.  Fig.  10  shows  a 

20 

few  of  the  low-lying  levels. 

The  first  experimental  observation  to  be  considered  is  the  decrease 
in  gas  temperature  and  pressure  immediately  following  the  input  of  energy 
from  the  laser.  This  seems  to  violate  laws  of  thermodynamics.  To  under¬ 
stand  this  effect,  examine  the  result  of  competing  energy  transfer 
processes  which  take  place  in  the  gas.  Ignore,  for  now,  effects  of 
thermal  conduction  to  walls  and  spontaneous  emission. 

The  laser  excites  the  mode  which  becomes  over  populated.  Absorp¬ 
tion  from  the  low-lying  level  to  level  can  also  occur  with  the 

same  result;  an  Increase  in  population.  Likewise,  a  molecules  in 
the  level  can  be  excited  to  2v^  +  v^.  From  Fig.  10  it  can  be 

seen  that  the  easiest  way  for  the  excited  molecules  at  level  to  reach 
equilibrium  is  to  undergo  a  v-v  transfer  to  an  energy  level  slightly 
higher  in  energy  (denoted  by  reaction  1).  The  energy  deficit  in  this 
reaction  must  come  from  somewhere.  It  comes  from  translation;  the  gas 
cools.  As  time  progresses,  energy  in  the  upper  vibrational  levels 
cascades  downwards  returning  energy  to  translation  resulting  in  heating. 
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Figure  10.  Partial  Vibrational  Energy  Level  Diagram  for  SF^. 
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Now  consider  what  will  happen  when  the  pressure  is  so  low  that 
thermal  conduction  and  spontaneous  emission  are  faster  than  the  transfer 
of  energy  to  translation  but  slower  than  the  Initial  v-v  process  which 
causes  the  gas  to  cool.  The  result,  shown  in  Fig.  8,  Is  that  the  gas 
initially  cools  but  rather  than  transferring  excess  energy  to  transla¬ 
tion,  it  is  transferred  to  the  tube  walls.  Since  the  test  cell  has  a 
large  specific  heat,  one  observes  only  cooling. 

Figure  11  shows  computed  curves  for  a  sequence  of  pressures  (0.30, 
0.25,  0.20  and  0.10  torr).  From  this  figure  it  can  be  seen  that  as  the 
pressure  is  decreased  the  amplitude  decreases  until  the  pressure  reaches 
0.2  torr  where  the  amplitude  changes  sign  (becomes  negative)  and  starts 
Increasing.  This  behavior  agrees  with  the  experimental  observation  shown 
in  Figure  8. 

Referring  back  to  Figure  7 ,  the  second  experimental  observation  is 
that  in  the  intermediate  pressure  regime,  the  curve  does  not  return  to 
zero  after  decaying;  it  overshoots  below  equilibrium.  This  Is  due  to 
acoustic  effects  which  have  not  yet  been  considered.  Figure  12  shows  the 
two  curves  computed  for  the  same  conditions  as  in  Figure  6  after 
Including  acoustic  effects.  It  can  be  seen  from  this  figure  that  the 
acoustic  effects  cause  the  curve  to  go  below  equilibrium,  returning 
slowly  to  equilibrium  on  a  longer  time  scale. 

It  can  be  seen  from  Figure  7,  also,  that  the  time  from  the  peak  of 
the  curve  to  the  point  where  it  goes  below  zero  is  0.40  msec  for  a 
pressure  of  0.4  torr  and  0.48  msec  for  a  pressure  of  0.80  torr.  Another 
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feature  of  this  figure  is  that  the  ratio  between  the  negative  part  of  the 
curve  to  the  positive  part  is  1/5  and  1/4  for  pressures  of  0.8  and  0.4 
torr,  respectively.  By  assuming  an  effective  radius  of  0.06  cm,  the 
computed  curves  in  Figure  13  fit  experimental  data  quite  well.  Figure  14 
compares  experimental  data  with  computed  points  for  the  low  pressure 
regime  (0.1  torr).  From  Figures  13  and  14,  it  can  be  seen  that  the 
computed  curves  agree  well  with  the  experimental  observations. 

Figure  15  shows  two  computed  curves  at  low  pressure  (0.05  torr).  The 
dashed  curve  shows  the  pressure  wave  when  spontaneous  emission  is 
included.  The  other  curve  shows  the  pressure  wave  when  spontaneous 
emission  is  excluded.  It  can  be  seen  from  this  figure  that  the  dashed 
curve  has  the  more  rapid  decay  time;  spontaneous  emission  offers  a 
parallel  pathway  for  deexcitation  processes  at  low  pressures. 

From  the  foregoing,  we  concluded  that  the  experimental  measurements 

can  be  reproduced  theoretically  by  assuming  a  relaxation  scheme  which 

connects  the  various  modes  of  SF,  collisionally  to  translation  with  rate 

o 

constants  for  the  various  processes  taken  from  other  measurements  or 
guessed.  The  rate  constants  were  varied  to  best  fit  the  data  in  Figure  7 
and  8  with  the  constraint  that  ultrasonic  and  laser  fluorescence 
relaxation  times  for  this  same  set  of  rates  agree  with  independent 
experimental  measurements  of  those  quantities.  This  same  procedure  was 
repeated  for  other  energy  transfer  pathways  until  best  agreement  between 
computed  and  measured  pressure  waves  was  obtained.  The  final  reaction 
scheme  assumed  is  given  In  Figure  10. 
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Figure  13.  Ftttlng  of  Experimental  Data  with  Computed  Points.  - >  Experimental 

Computed  Points. 
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Fitting  of  Experimental  Data  with  Computed  Curve  (Low  Pressure  Regime) 
Experimental  Data;  •••••,  Computed  Points. 


OS  torr 
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Figure  15.  The  Effect  of  Spontaneous  Emission  at  Low  Pressure.  - ,  with  spontaneous  emission 

••••,  without  spontaneous  emission. 
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